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Abstract 

Ideal hydrodynamic simulations are performed to compute the evolution with collision energy of 
hadron spectra and elliptic flow between AGS and LHC energies. We argue that viscous effects 
should decrease with increasing energy, improving the applicability of ideal fluid dynamics at 
higher energies. We show that the increasing radial flow at higher energies pushes the elliptic 
flow to larger transverse momenta, leading to a peaking and subsequent decrease of the elliptic 
flow at fixed pr with increasing collision energy, independent of whether or not there is a phase 
transition in the equation of state. 



1. Introduction and summary 

Ideal hydrodynamic simulations of the expansion stage of the hot and dense fireballs created 
in relativistic heavy-ion collisions predict a non-monotonic collision energy dependence of the 
(pr-integrated) elliptic flow V2(pr) fUl- The softening of the equation of state (EOS) at the 
quark-hadron phase transition leads to a predicted reduction of vi at RHIC energies, down from 
SPS energies, followed by another increase towards LHC energies. This effect is not seen in 
experiment [2] which shows instead a monotonic increase of with y[s. This is now understood 
as a failure of the ideal fluid picture during the late hadron gas stage which is highly viscous and 
inhibits the buildup of elliptic flow 13l|4]. Both viscous hydrodynamics 150 and hydro+cascade 
hybrid algorithms reproduce qualitatively the experimentally observed monotonic beam 

energy dependence of the integrated elliptic flow. Viscosity, in particular its strong increase in 
the hadronic phase, thus washes out the phase transition transition signature in the integrated 
elliptic flow excitation function. 

The PHENIX Collaboration observed that the pj -differential elliptic flow V2(pr), on the other 
hand, when plotted at fixed pr as a function of -\fs, shows signs of saturation at RHIC energies 
[6]. This has been interpreted as a possible remnant of the non-monotonic energy dependence 
predicted by hydrodynamics, signalling the softening of the EOS near T c and, possibly, even the 
existence of a critical end point (CEP) in the QCD phase diagram ]7[]. The apparent contradiction 
between a monotonically rising vii V*) as observed by NA49 and STAR yfl and a saturation 
withn increasing -\fs of the differential elliptic flow V2(pr) is resolved by the observation that also 
the radial flow increases monotonically with -\fs, leading to flatter /^-spectra at higher energies 
and thus pushing the hydrodynamically generated momentum anisotropy, which is reflected in 
viifr), to larger transverse momenta. 

Within a hydrodynamic picture of the collision fireball's collective evolution, the monotonic 
increase with -\fs of radial flow is a simple and unavoidable consequence of energy conservation, 
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independent of (and at most tempered by [ 1]) the existence of a phase transition in the QCD phase 
diagram. A systematic analysis of the hadron /^-spectra and viipr) as functions of collision 
energy [8] shows that a non-monotonic Vs-dependence of the elliptic flow V2(pr) at fixed pj, 
first rising from AGS to low SPS energies but then falling again towards RHIC and the LHC, is a 
generic consequence of the evolution of radial flow and, as such, cannot be used unambiguously 
as evidence in support or against the existence of the quark-hadron phase transition. To make 
this point is the purpose of this contribution. When searching for a clear QCD phase transition 
signature (in particular for the CEP), one has to look elsewhere. 

2. Ideal fluid dynamics from RHIC to LHC 

The analysis of Ref. [8] which is reported here is based on ideal relativistic fluid dynamics 
(IRFD). As already discussed above, IRFD is not perfect at RHIC energies and becomes increas- 
ingly worse at lower energies, due to the growing dynamical role played by the highly viscous 
hadron gas stage. At higher energies, the role of the hadronic phase decreases since more and 
more of the finally observed collective flow (in particular its anisotropy in non-central collisions) 
is generated already during the quark-gluon plasma (QGP) stage. The specific shear viscosity 
rj/s of the QGP (where s is its entropy density) is known to be very small, of the order of at 
most a few times the KSS [9] bound rj/s — 1 /An ifioll . For fixed t]/s, viscous effects in heavy-ion 
collisions are largest at early times, due to the large initial expansion rate from approximately 
boost-invariant longitudinal expansion. At any given early time t (before the onset of significant 
transverse expansion, t <k R/c s , where R is the transverse fireball radius and c s is the sound 
speed), viscous effects are controlled by the ratio of times scales -*• = - j-, where T s = rj/(sT) 
is the sound attenuation length and 1/r is the longitudinal expansion rate 11 111 . In perturbative 
QCD, the dimensionless specific shear viscosity rj/s is expected to increase only logarithmically 
with T Hence, r s is expected to decrease, leading (at the same r) to smaller viscous effects 
on hydrodynamic flow. Correspondingly, the validity of the IRFD approach should improve from 
RHIC to LHC. 

To extrapolate from lower to higher collision energies, we assume that thermalization occurs 
earlier at higher densities, i.e. at constant product 7oto = const. Using entropy conservation in 
IRFD, we can relate the final charged multiplicity to To and To as follows: dN c h/dy ~ dS jdy = 
to J d 2 x ± s(x±, to) ~ soto ~ toTq where so ~ is the peak value of the entropy density at To in 
central collisions. Combining both conditions we see that, starting from well-established initial 
conditions for 200AGeV Au+Au collisions at RHIC J3l, the initial thermalization time to and 

peak entropy density s scale as t ~ (^r)~ V2 , s ~ {^f-f 2 - The value of 45k for Pb+Pb at 
LHC energies cannot be predicted by hydrodynamics, but will be measured on the first day of 
LHC Pb-beam operation. We therefore present our results as a function of or, equivalently, 

of s - In H, the range s Q < 270 far 3 (^ < 1200) was explored; central 200AGeV Au+Au 

collisions at RHIC correspond to s = 1 17 far 3 and ^f- = 685. 

3. Results 

The left panel of Fig. [TJ shows the /?j--spectra of thermally emitted pions and protons (reso- 
nance decay contributions not included) as they evolve from low AGS to LHC energies. The flat- 
tening effects of increasing radial flow are clearly visible, especially for protons where strength- 
ening radial flow leads to a yield reduction at low pj in spite of the increasing total proton multi- 
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Figure 1 : (Color online) Left: Evolution of pion (top) and proton (bottom) transverse momentum spectra for central 
Au+Au collisions from low AGS to LHC energies (to correlate sq values with collision energies and charged hadron 
multiplicities, see Fig. 1 in Ql). Middle: Evolution of pion (top) and proton (bottom) differential elliptic flow in non- 
central Au+Au collisions at b = 7 fm. Right: Pion elliptic flow as function of collision energy (see text for discussion). 



plicity. The consequences of this shape change in the spectra for the /^-differential elliptic flow 
is seen in the middle panel of Fig. Q] again for thermally emitted pions (top) and protons (bottom) 
only. Initially, the differential elliptic flow V2(pr) increases from s = 8.3 far 3 (yfs ~ 4 GeV) 
to s = 23.4 far 3 (V* ~ 10 GeV), due the increase in total fireball lifetime before freeze-out 
which allows more elliptic flow to develop. At higher energies, however, increasing radial flow 
pushes the \'2{pt) curves to the right (more so for protons than for the lighter pions), leading to a 
decrease of elliptic flow at fixed pj. The bottom part of the right panel of Fig. Q]shows that (for 

> 10 GeV) this decrease of v^ 1 "" 1 PT {sfs) is monotonic, and that it holds for all pj values in the 
range pj < 1 GeV. Plotted logarithmically, the slope of this decrease is steeper for protons than 
for pions (not shown in Fig.[T|), reflecting the stronger radial flow effects on the heavier protons. 

The radial flow induced decrease of v^ sA PT (yfs) is independent of the behavior of the pj- 
integrated elliptic flow, shown in the upper part of the right panel of Fig. Q]for thermally emitted 
(dashed) and all pions (including resonance decays, solid). The integrated elliptic flow shows the 
well-documented non-mono tonic behavior of IRFD [1], featuring a decrease between top AGS 
and RHIC energies caused by the softening EOS near the quark-hadron phase transition, followed 
by an increase above -\fs > 100 GeV caused by the stiffening of the EOS in the QGP phase. The 
bottom panel shows that, at fixed pj, the differential elliptic flow continues to decrease while the 
integrated v'2 increases; these tendencies persist to the highest values of y[s where it is known that 
the elliptic flow fully saturates in the QGP phase, and that its finally observed value is therefore 
insensitive to the QCD phase transition and to the details of the conversion of quarks and gluons 
to hadrons. In this V^-region, it is obvious that the decrease of V2 Xed Pr (V*) is unrelated to the 
softening of the EOS near T c , and has therefore nothing at all to do with the phase transition. 

4. Conclusions 

Energy conservation and hydrodynamic behavior during the fireball expansion stage lead to 
increased radial flow from RHIC to LHC and correspondingly to flatter pj- and m^-spectra, 
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especially for heavy hadron species. As shown in Ref. JUl this causes baryon/meson ratios to 
continue to increase with both pj and mj-mo at LHC energies, as they do at RHIC. The slope 
of this increase as a function of transverse kinetic energy mj-rriQ is almost the same at LHC and 
RHIC, but as a function of pj the baryon/meson ratios increase with smaller slope at LHC than 
at RHIC, due to overall flatter /^-spectra. 

In ideal relativistic fluid dynamics (IRFD), the /^-integrated elliptic flow of pions and charged 
hadrons increases about 10-15% from RHIC to LHC energies; accounting additionally for vis- 
cous effects at RHIC (mostly of hadronic origin) that weaken or disappear at the LHC, the cor- 
responding increase is about 25%. At the same time, the differential elliptic flow at fixed pj, 
^fixed pt^^^ decreases from RHIC to LHC. This decrease is driven by an increase in radial flow 
which pushes the momentum anisotropy to larger pj\ it does not depend on a phase transition 
in the EOS. Combined with the increase of v 2 xe PT (yfs) at low -\fs < lOGeV, this leads to a 
non-monotonic -^-dependence G f y!| xed PT that is generic, caused by the interplay between radial 
flow and freeze-out, and not unambiguously associated with a phase transition in the QCD EOS. 
Although the analysis presented here was based on IRFD, the interplay between radial flow and 
freeze-out is a general principle that controls the buildup of elliptic flow also in real fluids. The 
observed non-monotonic energy dependence of PT is therefore robust, and (like variations of 
the EOS) inclusion of viscous effects is expected to only change the energy where PT peaks, 
but not the fact that it peaks. The search for QCD phase transition signatures, in particular for the 
predicted critical end point connecting a first order transition at high baryon density to a smooth 
cross-over transition at RHIC, cannot be based on this non-monotonic energy dependence of 
fixed-pr elliptic flow. 
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